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Abstract 

We analyze asymptotic scaling properties of a model class of anomalous reaction- 
diffusion (ARD) equations. Numerical experiments show that solutions to these have, 
for large t, well defined scaling properties. We suggest a general framework to ana- 
lyze asymptotic symmetry properties; this provides an analytical explanation of the 
observed asymptotic scaling properties for the considered ARD equations. 

Introduction 

In this note we will consider a class |15| of scalar partial differential equations (PDEs) of 
reaction-diffusion type associated to anomalous diffusion (see e.g. [3] for a recent review 
focusing on aspects of interest here). This is far from representing the most general 
anomalous reaction-diffusion (ARD) type of equation, but display a variety of behaviors 
common to much more general ARD equations. 

Numerical experiments on representatives of this class ^5] show that for large t solu- 
tions are described by travelling fronts with a well defined scaling behavior (see below for 
details); our goal is to provide an analytical explanation for this. 

In order to do this, we will first recall standard notions from symmetry analysis of dif- 
ferential equations (sect.l), and then extend them to the asymptotic framework (sect. 2). 
We will then be able to propose a general approach to extract asymptotic behavior of equa- 
tions based on maps to equations with known asymptotic symmetry properties (sect. 3); 
the basic idea will be, in the renormalization group langauge, to identify an equation in 
the same universality class amenable to asymptotic analysis. 

This approach will be used to analyze the asymptotic behavior of our class of anomalous 
reaction-diffusion equations in terms of the known asymptotic behavior of the FKPP 
equation. Our results provide a sound theoretical explanation of the behavior observed in 
numerical experiments ^3] and recalled below. 
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the work of RM was supported by INFM. - Istituto Nazionale di Fisica della Materia. 
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1 Symmetries of differential equations 

We assume the reader to be familiar with the main concepts and definitions for symmetries 
of differential equations (see e.g. ^3 El EE1 EE] for general treatments of these), to be 
later extended to asymptotic symmetries. In this section we will fix some general notation 
to be freely used later on, and recall general results we need later on; as in this note we 
will only consider scalar PDEs, we will specialize formulas to this case. 



General notation 

We consider an evolution PDE of order n for a real dependent variable u = u(x,t), with 
x, t independent real variables. We denote by M = R 2 x R the total space of independent 
and dependent variables, and by j( n >M the jet space of order n over M. Given a function 
u = f(x,t), we denote its graph as 7/ C M, and its prolongation as 7! C J^M. 
Let us consider a vector field 

d d d 

X = £{x,t,u) — + T(x,t,u)— + ip(x,t,u)— (1.1) 

in M. For e sufficiently small, the function u = f(x,t) is mapped by exp(eX) into a new 
function u = f(x,t); one obtains with standard computations that this new function is 
given by 

. c: f ) c en 



f = f + e 



dx J \ dt 



+ o(e) . (1.2) 

u=f(x,t) 



The action of X in j( n >M is described by its prolongation X^ n \ We say that X is a 
symmetry of a given equation if it maps any solution into a (generally, different) solution 
[HliniiniEl]- If the equation is given by $ = 0, this is equivalent to [X^(^)]^ =0 = 0. 

Geometrically, A identifies a solution manifold S in j( n 'M; the function u = f(x,t) 
is a solution to A if and only if 7^ C S, and X is a symmetry of A if and only if 
: S -» TS. 



Invariant solutions 

If A, which we write as <3?(x, t, u, ...) = 0, is a PDE for u = u(x, t) admitting a vector field 
X as symmetry, there is a well known procedure to determine X-invariant solutions to A; 
this represents an extension of the familiar method of characteristics to solve quasilinear 
PDEs (for more details, see e.g. the discussion in chap. 3 of [T7]). 1 

First of all we pass to symmetry-adapted coordinates (y, v, a) in M; the (y, v) will be 
X-invariant coordinates, while a will be acted upon by X. We see (y, a) as independent 
variables and v as the dependent one; we can then use the chain rule to express x and 
t derivatives of u in terms of the a and y derivatives of v, and write A in terms of the 
(a, y, v) coordinates, i.e. in the form $(0", y, v, ...) =0. As X is a symmetry of A, it follows 
that the function $, when subject to the side condition dv/da = 0, is independent of a. 

If we are able to determine a solution v = f(y) to the reduced equation, by writing this 
in terms of the (x, t, u) coordinates we get a X-invariant solution to A. 

1 This method requires a transversality condition, generically satisfied for scalar equations as those we 
wish to consider; see [TllLS) for discussion and a more general approach. 
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Equivalent equations 

Let us consider a map x '■ M — > M which is not a symmetry of A := Ao; we write 
X ■ (x, t, u) i — > (y, s, w) and x '■ A i— > A. 

If A is of the form A = ut — F(x,t,u,u x ,u xx ) = and x is projectable i.e. such 
that s = s(t), y = y(x,t), then A is of the same type: A = w s — G(y, s,w,w y ,w yy ) = 0. 
Moreover, if x is projectable it maps 7/ C M into a manifold which also is the graph of a 

function w = g{y, s), say x{lf) = 7g C M; this extends to prolongations, i.e. x^ : if^ ~~ > 

(2) ^ 
7g . Hence solutions u = f(x,t) to A are mapped into solutions w = g(y,s) to A. We 

say therefore that x is a solution preserving map |2U| . 

If x is invertible (with a C n inverse), we can repeat these considerations for x _1 - in this 

case the equations A and A of order n are equivalent, in that there is a C n isomorphism 

between solutions to A and solutions to A. 

Conditional and partial symmetries 

As mentioned above, A is a symmetry if it maps any solution to a (generally, different) 
solution. However, there can be cases where this is true only for some solution; one 
can formulate correspondingly weaker notions of symmetry. (For a review of different 
"extended" notions of symmetry, see fll]; see also for a recent and shorter discussion.) 

We say that A is a partial symmetry for A if there is a nonempty set Sx of solutions 
to A which is mapped to itself by A. The set Sx could be made of a single solution, and 
more generally that there could be solutions which are left invariant by X. In this latter 
case, we say that A is a conditional symmetry for A. 

It follows from (1.2) that u = f(x,t) is invariant under A if and only if 

A x := tp[x,t,f(x,t)]-u x £[x,t,f(x,t)]-UtT[x,t,f(x,i)] = 0; (1.3) 

thus A-invariant solutions to A are solutions to the system made of A and Ax- In other 
words, A is a conditional symmetry of A if and only if it is a symmetry of this system 
(and is not a proper symmetry of A). 

Partial symmetries can be seen in a similar way. As discussed in p[], a function / in 
the globally invariant set of solutions to A : $ = (/ 6 Sx in the notation used above) 
will be a solution to a system 

$(°) = , = , , = , 

where $(°) = $, and ¥ k+ ^ ■= The inte ger p, i.e. the order of the system, is 

determined as the lowest order such that $^ vanishes identically on common solutions 
to all the previous equations. Each equation <I>( fc ) = can, and should, be simplified by 
taking into account the previous equations; for concrete examples, see jjj. 

If A and A are equivalent via a solution-preserving map Xi this entails a corresponding 
relation between their conditional and/or partial (as well as ordinary) symmetries. 

2 Asymptotic symmetries of PDEs 

The concept of symmetry of a differential equation can be extended to that of asymptotic 
symmetry; we will again confine ourselves to scalar second order evolution PDEs. Asymp- 
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totic symmetries (in the present sense) were introduced and discussed in quite a general 
framework in ^U|, see also [5]. 

As suggested by their names, asymptotic symmetries of a differential equation A are 
vector fields X which, albeit in general not symmetries of A, satisfy X^ : Sa TSa 
asymptotically (see below for the precise sense of this). Any exact symmetry is also 
(trivially) an asymptotic symmetry. 

We write the equation A in the form 

$ = u t - F(x,t,u,u x ,u xx ) = , (2.1) 

and denote by T the space of maps F : (x,t,u,u x ,u xx ) — ► R which are polynomial in 
(u,u x ,u xx ) and rational in (x,t). Note that this space, which corresponds to the space of 
equations in the form we are considering, is invariant under scaling transformations and 
under translations. 

Vector fields will be written as in (1.1). The second prolongation of X will be denoted 
as as seen above, if X is projectable, then Y acts in T jllj . 

We write the flow (in the space of functions / : R 2 — ► R) issued by fo(x, t) as f\(x, t) = 
exp[AX], with df x /d\ = X[f\] = [tp - f x £ - f t T\ u =f( x ,t), see (1.2). We say that fo(x,t) is 
X-invariant if it is a fixed point for the flow of X, i.e. if X [/o] = 0. 

It may happen that /o is not X-invariant, but the flow fx is asymptotic to an invariant 
function i.e. 

lim \f x (x,t)-U(x,t)\=0 , X[f*} = 0. (2.2) 

A^oo 

When (2.2) is satisfied, we say that / is asymptotically X-invariant under the flow of X. 

Let us now consider the action of X, more precisely of Y, on the space of equations in 
the form (2.1), i.e. in T . We write 

A A = e Ay A := cr(A) [u t - F x (x,u,u x ,u xx )] (2.3) 

for the Y flow issued from A = Ao; by construction, this satisfies dA\/d\ = Y(A\). In 
this way, X induces a vector field W in the space T (W is nothing else than the restriction 
of Y to F); thus (2.3) is equivalent to dF/dX = W{F). 

As recalled above, X is a symmetry of A if and only if Y : S — > TS. This condition 
is now rephrased in terms of T by saying that X is a symmetry of Ao given by Ao := 
ut — Fq(x, t, u, ...) if and only if Fq is a fixed point for the flow of W. 

Suppose now that X is not a symmetry of Ao - i.e. Fq is not a fixed point for the flow 
of W in T - but that the flow F\ issued by Fq under e xw satisfies 

lim \F X -F*\ = , W(FJ = . 

In this case we say that X is an asymptotic symmetry for Fq, i.e. for the equation Ao- 
By construction and by the invertibility of the map exp(AX) - hence of exp(Al") - for 

A finite, if u = f(x,t) is a solution to Ao then u = f\(x,t) will be a solution to A\; and 

conversely if u = f\(x,t) is a solution to A^ then u = f(x,t) is a solution to Ao- 

In the limit A — » oo the invertibility of (2.3) fails if the flow goes to a limit point; we 

can nevertheless still say that solutions to the original equation flow into solutions to the 

asymptotic equation, provided both limits /* and A* exist. 

If the limit A* exists, it captures the behavior of A for large A, i.e. for large (or small, 

depending on the sign of a and b) t and \x\. 
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Solution-preserving maps and asymptotic symmetries 

The use of this construction, combining with solution-preserving maps, is the following. 

If there is an equation A whose asymptotic behavior is well understood, and such that 
there exists a solution-preserving map \ with %( 2 ) : A — ► A, we can study the asymptotic 
behavior of solutions u(x, t) to A by means of the asymptotic behavior of solutions w(y, s) 
to A. 

If x nas a smooth inverse, and W is mapped by x^ m ^o W, we can study the flow of 
W using 

e XW = [x (2)]-l Q e XW Q x (2) (2 4) 

Note that (2.4) remains valid also in the limit A — > oo, i.e. for the asymptotic regime. 

Needless to say, this approach is particularly convenient when the asymptotic behavior 
of the solutions w(y, s) to Ao - i.e. the behavior of solutions to A* - is simple and/or in 
some way universal (even better if Ao is a fixed point under W). 

Denoting as w*(y,s) the limit expression for the solutions to Ao, the asymptotic be- 
havior of the solution u(x, t) to Ao will be given by 

u*(x,t) = x' 1 [w*(y,s)} . (2.5) 

More precisely, considering jf and j g corresponding to u = f(x,t) and w = g(y,s) (see 
sect.l), and going asymptotically into and 7* respectively, we have 

7/ = X- 1 (7*) ■ (2.5') 

Conditional and partial asymptotic symmetries 

It is also possible to consider asymptotic versions of conditional and partial symmetries. 

Let the function u = fo(x,t) (more precisely, the corresponding graph 70 C M) flow 
under e xx to a fixed point u = f*(x,t) (more precisely, a A-invariant graph 7* C M) 
albeit Ao does not flow to a fixed point. 

In such a situation, the solution manifold S\ C J^M does not go to a limit manifold, 
but there is a submanifold C S\, with 7^ C C 5a, which flows to a fixed limit 
submanifold , with 7* C . In this case we say that A is a conditional asymptotic 
symmetry for A. 

The same construction, with suitable and rather obvious modifications, applies for what 
concerns partial symmetries. Suppose that A = A^ ) does not flow to a fixed point under 
W, and consider the equations 

A (i) := y[ A (°)] , ... , AW : = YfA^'- 1 )] 

up to an r - if it exists - such that A^ r ) does admit a fixed point under the W flow, 
while the A^ with k < r do not. Then the manifold 

sf n...n5f := s x 

(with C J^M the solution manifold for A^) flows to a limit submanifold 5*, and 
solutions u = fo(x,t) to the system A^ (k = 0, ...,r) flow to functions u = f*(x,t) such 
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that the prolongation 7* C J^M of the corresponding graph 7^ = {(x, t, f*(x, t))} lie 
in S*. In this case we say that X is a partial asymptotic symmetry for A. 

If A and A are related by a solution-preserving map x^ 2 \ then x also relates their 
conditional and partial symmetries. In particular, if \ is invertible and A admits X as a 
conditional symmetry, then A admits x* 1 ^) as a conditional symmetry, with x* 1 the 
extension of x^ 1 to vector fields. This will be of use below. 

3 Asymptotic symmetries as a tool to test asymptotic 
behavior 

In many physically relevant cases, one has to study nonlinear PDEs which are not amenable 
to an exact treatment, or at least for which such a treatment is not known, and for which 
a numerical study shows an asymptotic behavior which appears to be well described by 
some kind of invariance. Usually, the latter corresponds to a scale invariance (self-similar 
solutions), or a translation invariance (travelling waves), or a combination of both of these. 

The discussion conducted so far can be implemented into a procedure allowing on the 
one hand to test if the observed asymptotic behavior is a characteristic of the equation 
(rather than an artifact of the numerical experiments conducted on it), and on the other 
hand to formulate simpler equations extracting the asymptotic behavior. 

We will now describe the procedure in operational terms; as a (rather simple, but 
relevant) example to illustrate our procedure we will consider here the heat equation, and 
then anomalous diffusion equations, while in later sections we apply the procedure on 
anomalous reaction-diffusion equations. 

We denote by X the vector field describing the observed invariance, and consider a 
second order equation for u = u(x,t) of the form A := ut — F(x,t,u,u x ,u xx ) = (hence 
M = {(x,t;u)}). 

• Step 1. Pass to symmetry-adapted coordinates in M, i.e. coordinates (a,y;v) such 
that X{y) = X(v) = 0; thus in these coordinates X = f(a,y,v)(d/da). 2 

• Step 2. Identify v as the new dependent variable, i.e. v = v(a,y). This allows 
to write x and t derivatives of u as a and y derivatives of v, hence to write the 
differential equation A(x,t;u^) as A(cr, y; v^). 

• Step 3. Reduce the equation A(a,y;v^) to the space of X-invariant functions, i.e. 
to v(a,y) satisfying v a = 0. For X an exact symmetry, the reduced equation Ax 
will not depend on a at all. For X a conditional or partial symmetry, a will still 
appear parametrically in the reduced equation. 

• Step 4. Study the asymptotic behavior of the solutions to the reduced equation 
Ax for a — > 00. 

• Step 5. Go back to the original variables. 



2 We stress that we are not requiring / = 1; actually when we deal with scaling symmetries it is 
appropriate to require f(a, y, v) = a. 
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Elementary example: the heat equation 



Let us illustrate our procedure by applying it on the heat equation ut = u 
solutions are of the form 



u(x, t) 



A 



cxp 



Ax 



21 



and are invariant under the scaling vector field 



X 



vd x + 2tdt — ud u 



Its asymptotic 
(3.1) 

(3.2) 



which is also an exact symmetry of the equation. 

The symmetry-adapted coordinates are a = t, y = x 2 /t, and v = xu. In these, the heat 
equation reads 

av a = 4:yv yy + (2 + y)v y + v. (3.3) 

Note that (3.3) necessarily requires that for a — > oo, v a = 0: this shows in a simple way 
that the full asymptotic behavior of (3.3) is captured by (3.1). 

Imposing v a = in (3.3), a disappears completely. Needless to say, the equation 
obtained in this way has solutions v(y) = A^/y exp[— y/4], which when mapped back to 
the original coordinates produce the gaussian (3.1). 



Example: anomalous diffusion equations 

The procedure described above can also be applied to what will be our model class of 
anomalous diffusion equations, studied numerically in |15j : these are written as 



x- *l 2 d r i 



l-al- ,. 

x' Q /V 



t x - va dx 



(3.4) 



To focus ideas, let us mention two examples of equations in this class: (i) For a = 2 
we have the generalized gaussian process; (ii) For a = 2/3, v = 3/2 we deal with the 
Richardson equation describing the evolution of the distance between two particles in 
developed turbulent regime; we refer to ^3] for a discussion of the interest of the class of 
anomalous RD equations (3.4). 
One can check that the map 

s = t ™ , y = x a ' 2 , w = t( 2 - a ^^u (3.5) 

is solution preserving from (3.4) to the heat equation w s = w yy . Using the inverse 
change of coordinates, the universal asymptotic solution w(y, s) ~ s~ x l 2 e~ y ' I s of the heat 
equation is mapped back into 

u(x,t) ~ 1 e W [-x a /t»] , (3.6) 



which represents therefore the universal asymptotic solution to (3.4). This result is con- 
firmed by numerical experiments |15j . 
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4 The FKPP equation: asymptotic solutions, and 
symmetries 

In the same way as in the previous example the heat equation played the role of target 
for solution preserving maps applied to anomalous diffusion equations, in the case of (our 
model class of) anomalous reaction-diffusion (ARD) equations we will look for a solution- 
preserving map to the well known Fisher-Kolmogorov-Petrovskii-Piskunov (FKPP) equa- 
tion [HI EES EE]. This reads 

u t = Du xx + eu(l - u) , 

with e and D real positive constants, and one requires that u(x,t) > for all x and 
t. There are two stationary homogeneous states, i.e. u = and u = 1; the latter is 
stable while the former is unstable against small perturbations. In this section we study 
asymptotic symmetries of the equation and of its solutions. 

It is well known j(.4j that if the initial datum is suitably concentrated, e.g. u(x,0) = 
for \x — xo\ > L or more generally u(x,0) < Aexp[—x/L], then asymptotically for t — » oo 
and x — > oo the solution is of the form u = f(x,t) ~ exp[— (x — vt)/X], with A = \/D/e 
and v = V AeD. This represents a front of width A travelling with speed v; it connects the 
stable state u = 1 and the unstable state u = 0. 

In discussing the FKPP equation, it is convenient to pass to rescaled coordinates t = et, 
x = (y e/D)x. From now on we will use these coordinates, and omit the tildas for ease of 
notation. In these coordinates, the FKPP equation reads 

Ut = u xx + u(l - u) . (4.1) 

As for the asymptotic solution described above, this reads now 

u = f (x,t) ~ A exp[-(x-2t)] ; (4.2) 

note the front has speed v = 2 and width A = 1. 

It should be stressed that the f(x,t) or fo(x,t) given above provide the solution for 
x — ► oo, i.e. in the region of small u; in this region (4.1) is well approximated by its 
linearization around u = 0, i.e. 

u t = u xx + u ; (4.3) 
the ansatz u(x,t) = w(z) := w(x — 2t) takes this into the ODE 

w" + 2w' + w = (4.4) 

for w = w(z), with solution 

w(z) = cie~ z + c 2 ze~ z . (4.5) 

We denote by W the set of solutions described by (4.5); note that W = R 2 , and (01,02) 
provide coordinates in W. 

The fo given above, see (4.2), corresponds to C2 = 0. This can be characterized in 
terms of symmetry properties, as discussed below. It is convenient to consider linear 
combinations of the shifts, given by X± = X\ =p (1/2)^2 ; note that z = x — 2t is invariant 
under X_, and that X + = d z . We also have Xq = wd w . Needless to say, [Xq,X + ] = 0. 
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Lemma 1. The symmetry algebra of the linearized equation (4-3) is generated by the scal- 
ing Xq = ud u , and by the translations X\ = d x and X2 = dt- The quotient equation (4-4) 
admits only Xq as scaling symmetry; it also admits the translation symmetry generated by 
X + , while X- has been quotiented out by passing to the z variable. 

Proof. This follows from standard (and elementary) computations. 

Lemma 2. The propagating front solutions correspond to a subspace of'W invariant under 
the action of the group generated by the vector fields Xq and X + . 

Proof. A general element of the group generated by Xq and X + is written as g(a,[3) := 
exp[aA~o + (3X + ] and acts on W by 

g(a,/3) : ( Cl ,c 2 ) - (e a+ ^( Cl + (3c 2 ) , e Q+ ^c 2 ) . 

The subspace c 2 = is invariant under this action. As remarked above, the propagating 
front solutions (7.2) correspond to c 2 = 0. 

It is immediate to see that the only scaling or shift symmetries of the full FKPP equation 
(4.1) are those, with generators X\ = d x and X2 = dt, corresponding to translations in x 
and t; these reflect the fact that (4.1) is a homogeneous equation. 

The situation is different for what concerns asymptotic symmetries, and in particular 
scaling ones, as we now discuss. 

Lemma 3. Let X be a scaling vector field, such that lim^oo exp(AA) extracts the behavior 
for large \x\ and t. Let Aq be the FKPP equation, and A\ = e XY Aq with Y the prolongation 
of X. Then lim^oo A\ = A* is the heat equation ut — u x 



Proof. We consider the most general scaling generator, i.e. a vector field in the form (1.3), 
X = axd x + btdt + cud u . We can always set one of the constants (a, b, c) equal to unity 
(provided it is nonzero); this amounts to a redefinition of the scaling group parameter. 

Applying the procedure described in previous sections, with of course Ao := ut — u xx — 
u(l — u) = the FKPP equation, we obtain at once that 



u t 



\ b ~ 2a u xx - \ b u + X b+C u 2 



(4.6) 



We choose c = b and a = b/2. In order for lim^oo exp(XX) to extract the behavior for 
large \x\ and t, we must choose a < 0, b < 0. We can set the modulus of one of the 
constants, say b for definiteness, equal to unity; i.e. b = — 1. With these choices, we have 

A x = in - u xx - A _1 u + A"V . (4.7) 

The limit A* := lim^oo A\ is the heat equation ut — u xx = 0, as claimed. 



5 Anomalous reaction-diffusion equations 

The FKPP is a (relevant) representative of a more general class of anomalous reaction- 
diffusion equations, which we write as 

u t = L[u] + h(u) ; (5.1) 
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here L is the linear operator describing passive transport of the field u, hence anomalous 
diffusion, while h(u) describes its growth. Logistic growth, which we will consider, corre- 
sponds to choosing h(u) := u(l — u). With this choice and considering the L associated 
to anomalous diffusion in our model class, see the r.h.s. of (3.4), we get the equation 



ut 



„2-a 



t l-ua 



(2 -a) ' 

Uxx ~\~ 2 ^ x 



+ u(l-u) . 



(5.2) 



One is usually interested in solutions with initial data u(x, 0) which are suitably regular 
and with compact support. A detailed numerical study of systems of the form (5.2) with 
such initial data was conducted in |15j . 

We summarize the findings of these numerical experiments as follows: 

• (i) asymptotically for large x and t, the solution is described by a travelling front 
with varying speed c(t) and width A(i); the form of this front for small u is well 
described by 

x — c(t) ■ t 



u(x, t) ~ A exp 



X(t) 



(ii) the (asymptotic) time dependence of c(t) and X(t) are described by 



c(t) ~ c • t s , \(t) 
where Co and Ao are dimensional constants; 
(iii) the scaling exponent 8 is given by 

5 := v + (1/a) 



(5.3) 



(5.4) 



Thus we rewrite (5.3) in the form 

u(x, t) ~ A exp 



x - (cot^t 1 
A^F" 



(5.5) 



(5.6) 



Note that for 5 = 0, i.e. for v = 1 — (1/a), the front travels with constant speed and 
width, as for the FKPP equation. 

We want now to describe precisely the invariance properties of the observed asymptotic 
solution (5.6), in particular for what concerns scaling transformations. 

Lemma 4. The scaling invariance of (5.6) is described by the generalized scaling group 



x -» (X s ) x 
with X the group parameter. 



£->• Xt, [exp ((A- l)Kt)\ u 



(5.7) 



Proof. The generator of the one-parameter group described by (5.7) is 

X = 5xd x + tdt — Ktud u . 



(5.8) 



The invariance of (5.6) under this can be easily checked using (1.2). It can be shown that 
this is the only scaling type symmetry of (5.6). <0* 
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Theorem 1. The vector field (5.8) is not a symmetry of the equation (5.2), but it is an 
asymptotic symmetry of the same equation. 

Proof. Let us denote (5.2) by Ao and (5.8) by A; applying Y = on Ao, and 

restricting to the solution manifold So of Ao (this amounts to substituting for w a according 
to Aq itself), we obtain 



Ai := [y(A )] 5n = [(1 - a)t aS (x/t) 2 ~ a } u xx + 



\2-on 



+ 



(l/2)(a - l)(a - 2)t aS ~ 1 {x/t) l - a 
u(l + K - u + Kut) . 



u x + 



This is not zero, i.e. X is not a symmetry of (5.2). 

Going further on with our procedure, we have to compute A2 := [y(Ai)] s nSi . It 
results 

A 2 := [2 - AKt + K 2 t 2 - a(l - Kt)} u 2 - [(1 + K){a - 2)] u = 0. 

This has the trivial solution u = 0, which is also solution to Ao and Ai, and the nontrivial 
solution 

u(t) = (2-a)(l + K) 

{ ' (2 - a) - (4 - a)Kt + K 2 t 2 ' 

The latter, as easily checked by explicit computation, is in general not a solution to Ao 
and Ai: inserting this into Ao and Ai we have respectively 

'Kt 2 + (2K - 4 + a)t + (2a - 6) ' 



A = [{a-2)K{l + K)\ 



(K 2 t 2 + (a- 4K)t + (2 - a)Y 
Kt 2 - 2t 



Ai = \{a-2)K(l + K) 2 } , 

J \{K 2 t 2 + (a - AK)t + (2 - a)) 2 

For K 7^ 0, —1, both of these expressions are not zero (unless a = 2, corresponding to 
gaussian processes). However, both of these go to zero (like 1/t 2 ), for all a and K, in the 
limit £ — > 00. 

Having determined that A is an asymptotic (partial) symmetry for our equation Ao, 
we will now apply our general procedure. The first step consists in passing to symmetry 
adapted coordinates; these are 



a = t , y = x/t , w 



ue 



(5.9') 



In these coordinates, the vector field (5.8) reads simply A = crd a , and the (obviously 
A-invariant) asymptotic solution (5.6) is w = Aexp[y/\o\. With standard computations 
we obtain 



U-r 



w a - 5{y/o)w y - Kw] e Ka , 
\l/a s )w y ] e~ Ka 
\l/a s ) 2 w yy ] 



-Ka 



(5.9") 



Using (5.9') and (5.9"), we express Aq in the new coordinates; this results to be 



,2-a 



o-x 



Wyy + 



a 



1-a 



at* 



+ a 



w y + (K + l)w — e 



-Ka w 2 



(5.10) 
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where we have defined \x = a(5 — v + 1/a) for ease of writing. 

The expression (5.10) holds for the general map (5.9); however we are specially inter- 
ested in the choice 5 = (y — 1 + 1/a), see (5.5). With this, we have fi = (2 — a), and finally 
(5.2) reads 



■»' n 



a 



2-a 



Wyy + 



2 - a \ fyV~ a , y 



w y + (K + l)w - e- Ka w 2 . (5.11) 



In the limit a — > oo, the last term disappears (faster than any power in a), and (5.11) 
reduce to a linear equation. 

Theorem 2. The equation (5.11) has no nontrivial X -invariant solutions, but admits 
nontrivial asymptotically X -invariant solutions. 

Proof. The X-invariant solutions to (5.11) are obtained by requiring that w a = 0; with 
this the equation reduces to 



(!) 



2-a 



Wyy + 



2-a \ (y 
2a 



a 



l-Q 



y 

+ a— 
a 



w y + (K + l)w — e 



-Ka w 2 



. (5.12) 



Note that a appears parametrically here, and (5.12) splits into the equations corresponding 
to the vanishing of coefficients of different powers of a (this is a general feature of partial 
or "weak" symmetries, see jSj). The only common solution to these is w = 0, which proves 
the first part of the statement. 

Let us go back to considering (5.11); in order to study its asymptotic behavior for 
a — ► oo, we disregard the term which is exponentially small for large a. The resulting 
linear equation for w = w(y), i.e. 



a J 



2-a 



Wyy + 



2-a 
2a 



l-a 



+ a 1 



a 



w y + (K+l)w = , 



yields as solution 



■»' 



(y) = Ci£[0,0] + (^2(K + l)a2(y/a)«) C 2 £[1/2,1] 



where C\,C2 are arbitrary constants, and 

IC[x,y] := F u [{K + l)a/a 2 + x; 1/2 + y; -a(y/a) 
with Fu = \F\ the Kummer confluent hypergeometric function 

r(6) 



F u [a;b;z] := iFi(a;b;z) 



T(b-a)T(a) J 



e zt t a - l {l-t) b - a - v dt . 



b-a-l 



The asymptotic solution could now be expressed in terms of the original variables using 
(5.9); this yields an explicit but involved and not specially illuminating expression. <3> 
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6 Other asymptotic partial symmetries of ARD equations 

The scaling symmetry (5.8) is not the only symmetry of the observed asymptotic solution 
(5.6) to (5.2). In this section we identify different symmetries and apply our approach on 
the basis of these. 

Lemma 5. The vector field X = ^d x + rdt + ipd u is a symmetry of the asymptotic solution 
(5.6) to the equation (5.2) if and only if it belongs to the two dimensional module (over 
smooth real functions C°°(R 3 , R) ofx,t,u) generated by 

X 1 = 8 X + {J-^j d u ■ X 2 = (A t 1+5 ) d t - (x5 + cot l+s ) d u . (6.1) 



Proof. This follows easily by using (1.2) and the explicit expression (5.6) of the asymptotic 
solution u = f*(x,t). Indeed, applying (1.2) we get 

(/-/) „„_r,_, 1 (xS + CQt^T-tt 



ip — A exp (cot — xt )/Aq 



£ 

and the result follows immediately 
In the following, we will consider in particular 

X := (x5 + c i 1+<5 ) d x + td t , (6.2) 

as well as X\ and X 2 themselves. Note that Xq operates in the space of independent 
variables alone. 

We will write second-prolonged vector fields in the form 

r\ r\ r\ r\ r\ 

Y = X< 2 > = X + V x — + + V xx - + + ^ttj— ; (6.3) 

OU x OUt <JU XX UU x t UUft 

in view of (5.2) we will actually need only the coefficients ^/ x , ^ xx . We will, as in the 
previous discussion, denote (5.2) as Ao- 

Theorem 3. The vector fields Xq, X\ and X 2 are partial symmetries of A . 

Proof. We will denote by Y{ the second prolongations of X{. Let us start by considering 
X\. In this case the coefficients of the second-prolonged vector field Y\ are: 

* x = , * t = -(<yAo)t~ (1+5) ) , = , V xy = 0, = (<5/A )(l + 5)t-^ . 

We define then Ai = [Yi(Ao)] s ; A2 = [Yi(Ai)] 5(jn5i . By explicit computation, it results 
that 

this shows that X\ is a partial symmetry for Ao- 

For the other vector fields, note that the relevant coefficients of Y 2 are 

5 (l + S)Sx T 
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and those of Yq are 

*$> x = -5u x , * t = -(1 + 5)cot 5 u x - u t , *$> xx = -25u xx . 

Using these, the theorem follows by explicit computations. 

It turns out that reduction, and invariant solutions, under the vector field Xq are 
of special interest. This is due to the following theorem, which provides an analytic 
explanation of the numerically observed behavior. 

Theorem 4. The equation (5.2) admits an asymptotically X^-invariant solution, described 
by (5.6). 

Proof. In this case the symmetry adapted coordinates are 

a = t , y = (x/t s ) — cot , w = u ; 
the relevant u derivatives are expressed in the new coordinates as 

f Sy + cq(1 + S)a\ 1 1 

U t = W a - I \Wy, U X = JgWy , U XX = ^Wyy . 

In these coordinates the equation (5.2) is written as 

W a = AWyy + BWy + f (w) , (6.4) 

with 

A = (V + W) 2 ~ a ; B = (co + S y -±^ + e( y -±^-X~ a ^— ) . 
\ a J y a \ a J y + c a J 

The vector field Xq reads simply Xq = ad a ; its second prolongation is 

d d d d 

Yq = a- W a - W ay - ^W ua - . 

oa ow a ow a y ow aa 

For o" — > oo, (6.4) reads 

W a = cl~ a Wyy + C0(l + 5)Wy + f (w) . 

The Xo-invariant solutions satisfy w a = 0, and are thus obtained as solution to 

cl~ a Wyy + C (l + 5)Wy + W = , (6.5) 

where we have used w « 1 in the region we are investigating (i.e. for a — > oo), so that 
f(w) ~ w. 

Solutions to (6.5) are of the form 

w(y) = Cie" u)+y + C2e~ u) ~ y , 

where 



1± wi 



c%(i + sy 
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If we require the solutions to be non oscillating, this implies a lower bound on the parameter 
c , i.e. co > (2/(1 + <5)) 2 /«. 

The solution e~ UJ+z is unstable against small perturbations, while e~ ul ^ z is stable |14| . 
As proved by Kolmogorov, the asymptotic solution is the stable one with the lowest speed 
giving nonoscillating behavior, i.e. Co = [2/(1 + 5)] 2 ^ a . This means w(y) ~ e~ u ° y with 
ujq = [2/(1 + 5)] 1-2 / a . Going back to the original variables, we get 



u(x,t) ~ A exp 



x — v(t)t 



A exp 



-LO - 



X — Cot 



1+8 1 



t s 



This is precisely the numerically observed asymptotic behavior (5.6). 



7 Conclusions and discussion 

We provided suitable definitions of asymptotic symmetries - both in proper and in condi- 
tional or partial sense - and proposed a method for the analysis of asymptotic symmetry 
properties of PDEs and their solutions. 

We applied our general method to a model class of anomalous reaction-diffusion (ARD) 
equations, discussed and studied numerically in |15j . We first considered the standard 
FKPP equation, and described in detail its asymptotic symmetry properties; we have 
also shown that our approach recovers the well known asymptotic properties of FKPP 
solutions. 

We have then tackled general ARD equations in our model class, i.e. with anomalous 
diffusion associated to (3.5). We recalled the features of asymptotic solutions as observed 
in numerical experiments, and identified the Lie generator X of the observed asymptotic 
generalized scaling invariance; in theorem 1 we showed that this is not a symmetry, but it 
is an asymptotic symmetry, of the ARD equation. We have then considered the solution- 
preserving maps associated to this asymptotic scaling symmetry, focusing on the physical 
value of the parameter 5; in theorem 2 we have shown that in this case the ARD equation 
has no solution invariant under X (which therefore is not a conditional symmetry of the 
equation), but admits solutions which are asymptotically invariant under it. That is, X is 
an asymptotic conditional symmetry of the ARD equation. 

Finally, in sect. 6 we passed to consider in more detail the numerically observed asymp- 
totic solution (5.6). We identified the full symmetry algebra Q of it; vector fields in this 
algebra are asymptotic conditional symmetries of the ARD equation. We focused in par- 
ticular on certain vector fields in Q, and shown they are partial symmetries for the ARD 
equation. Among these vector field is the scaling vector field Xq given by (6.2) and not 
depending, nor acting, on the dependent variable u. We proved that the ARD equation 
does admit an asymptotically Ao-invariant solution, which is precisely the numerically 
observed one (5.6). 

We have thus provided an analytic explanation for the numerically observed behavior, 
based on our general method. 

Let us conclude by presenting some brief final remarks on our method, 
(a) This method represents an evolution of the classical method to determine partially 
invariant solutions for symmetric PDEs |17[ 121 j , and a blend of it with the method of 
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conditional and partial symmetries [Zj, in order to analyze equations which do not have 
complete (as opposed to asymptotic) symmetries. Our method is based on the abstract 
approach developed in ^JJl, based itself on ideas and previous work by several authors 

iiiiaiEid]. 

(b) The application of the method to (generalized) scaling symmetries and our model 
class of ARD equations was greatly facilitated by the form of the vector fields and of the 
initial equations Aq: indeed, the ^-accessible part of T was finite-dimensional. 

(c) Our method deals, strictly speaking, with properties which are asymptotic in the 
group parameter; these correspond to properties asymptotic in time and space only if the 
considered vector field has favourable properties itself. Also, our method does not intend 
to tackle intermediate asymptotics 

(d) This work was concerned only with scale invariance (at infinity or near a travelling 
front). We trust however that suitable generalizations of our approach can also deal with 
more general asymptotic invariance properties, and more general differential equations, 
and that it is potentially capable to provide a sound explanation - or prediction - of the 
asymptotic invariance of their solutions. 
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